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Abstract 

Modern developments in population dynamics emphasize the role of the turnover 
of individuals. In the new approaches stable population size is a dynamic equilib- 
rium between different mortality and fecundity factors instead of an arbitrary fixed 
carrying capacity. The latest replicator dynamics models assume that regulation 
of the population size acts through feedback driven by density dependent juvenile 
mortality. Here, we consider a simplified model to extract the properties of this ap- 
proach. We show that at the stable population size, the structure of the frequency 
dependent evolutionary game emerges. Turnover of individuals induces a lottery 
mechanism where for each nest site released by a dead adult individual a single 
newborn is drawn from the pool of newborn candidates. This frequency dependent 
selection leads toward the strategy maximizing the number of newborns per adult 
death. However, multiple strategies can maximize this value. Among them, the 
strategy with the greatest mortality (which implies greater instantaneous growth 
rate) is selected. This result is important for the discussion about universal fitness 
measures and which parameters are maximized by natural selection. This is related 
to the fitness measures Rq and r, because the number of newborns per single dead 
individual equals lifetime production of newborn Rq in models without ageing. Our 
model suggests the existence of another fitness measure which is the combination 
of Rq and r. According to the nest site lottery mechanism, at stable population 
size, selection favours strategies with the greatest r from those with the greatest 

Rq. 

Keywords: density dependence, logistic equation, frequency dependent selection, 
life history, evolutionary game, fitness measures 
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1. Introduction 

In the current development of th evolutionary ecology (Post and Palkovacs 2009, 
Pelletier et al. 2009, Morris 2011, Schoener 2011) a problem of the eco-evolutionary 
feedback is of the special interest. One of the major theoretical problems in the 
modelling of population dynamics, and in general of evolutionary biology and ecol- 
ogy, is the limit of population growth and its selection consequences. This topic is 
very important in many disciplines such as evolutionary game theory and life history 
theory. The earliest attempt to solve this problem for populations with overlapping 
generations is the continuous logistic equation introduced by Verhulst in the 19th 
century (Verhulst, 1838), which can be found in every textbook on ecology and 
mathematical biology. It inspired the idea of r and K selection (McArthur and 
Wilson 1967), that selection favours different strategies at low densities and near 
the stable population size, and is still applied in modelling (Cressman et al. 2004, 
Cressman and Krivan 2006, Cressman and Krivan 2010). This concept states that 
there is some arbitrary maximal population size at which growth is suppressed and 
the population remains stable. However, this approach produces some unusual pre- 
dictions which provoked a wide discussion (Kozlowski 1980, Lomnicki 1988, Kuno 
1991, Ginzburg 1992, Gabriel 2005, Hui 2006, Argasinski and Kozlowski 2008) pre- 
sented in the next section. The problem of the limits to growth is important not 
only for ecological population growth models but also for the modelling of natural 
selection. There is huge discussion on what is maximized by natural selection and 
what happens when a population reaches the limit to growth (Metz et al. 1992, 
Kozlowski 1993, Mylius and Diekmann 1995, Brommer 2000, Dieckmann and Metz 
2006, Metz et al. 2008, Roff 2008). However in these attempts density dependence 
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is represented by some abstract unspecified factors. Thus tlie proposed solutions are 
very general and abstract. A concrete mechanistic interpretation should be help- 
ful in the interpretation of the mathematical notions. In the modification of the 
logistic equation (Kozlowski 1980, Hui 2006) which was applied in game-theoretic 
modelling (Argasinski and Kozlowski 2008, Argasinski and Broom 2012) there is an 
example of a mechanism responsible for strategically neutral density dependence 
called in this paper a "nest site lottery". The underlying assumption is that there 
is a limited number of nest sites at the environment and that every newborn must 
find a nest site to survive. Thus all newborns produced in some time moment 
form a pool of candidates to be drawn from to replace the dead individuals in their 
nest sites. The difference is that in this case there is arbitrary maximal population 
size described by a carrying capacity indicating the number of available nest sites 
(mechanistically interpreted as a nests or holes, where individual can settle, as in 
Hui 2006). However, the stable population size is not the carrying capacity, as in 
the classical logistic equation, but the dynamic equilibrium between different fac- 
tors of mortality and fecundity (Kozlowski 1980, Ginzburg 1992, Hui 2006), which 
can be affected by the dynamics of the population state (Argasinski and Broom 
2012). The advantage of this approach is that it considers a realistic turnover of 
individuals (Argasinski and Kozlowski 2008). In this paper we will more rigorously 
analyze the properties of the nest site lottery mechanism in a simpler model than in 
the previous papers (Argasinski and Kozlowski 2008, Argasinski and Broom 2012). 

Below we give the mathematical details of this approach (sections 2 and 3). Sec- 
tion 4 starts the development of the selection model and at section 5 the important 
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notion of the turnover coefficient will be introduced. Section 6 contains presenta- 
tion of selection dynamics. Section 7 contains the main result which are equations 
(|13I12I) and their analysis describing the nest site lottery mechanism ( intuitively 
depicted at Fig. 1). The mechanistic reasoning from section 7 is completed by 
Theorem 1 describing quantitative characterization of the restpoints of the system 
()13|12p . In section 8 the relationships between the nest site lottery mechanism and 
the invasion fitness concept are presented. The importance of the obtained results 
and general ideas inspired by them are presented in the discussion (the last subsec- 
tion contains discussion on the existence of the new fitness measure suggested by 
obtained results). 

2. Density dependence 

The cornerstone of mathematical ecology is the Malthusian equation describing 
exponential population growth, 

(1) — = nr = nb — nd = n (b — d) , 

dt 

where b is the number of newborns and d describes the fraction of dead indi- 
viduals per unit time. Thus the Malthusian coefficient r can be interpreted as the 
balance between mortality d and fertility b. The above model is not realistic, be- 
cause it allows for infinite population growth. The classical solution of this problem 
is the use of the logistic equation, which is equation ([2]), 



However, this relies on a problematic assumption which has very serious conse- 
quences. Equation ([2|) for example produces artifacts in population growth models 



(2) 
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(Kuno 1991, Gabriel 2005) and selection models related to replicator dynamics (Ar- 
gasinski and Kozlowski 2008). For example, it suppresses the selection dynamics in 
the replicator dynamics by setting the right hand sides of equations to (Argasin- 
ski and Kozlowski 2008), the trajectory escapes to infinity for r < (i.e. b < d) 
and initial size greater than K (known as Levins' paradox, (Gabriel 2005)) or the 
trajectory decreases with increasing rate for r < and initial size slightly smaller 
than K (Kuno 1991). This is caused by the fact that the whole term r is multiplied 
by the suppression coefficient, which implies that with population growth, both 
mortality and fertility decrease, and mortality decreasing with increasing popula- 
tion size and reaching zero at equilibrium is biologically counterintuitive. Mortality 
should not decrease with population growth and individuals cannot be immortal 
at equilibrium. The above problems suggest that models should rely on clear and 
mechanistic assumptions (Geritz and Kisdi 2012). Thus, density dependent sup- 
pression should act only on the number of juveniles recruited to the population 
(Kozlowski 1980, Ginzburg 1992) and the initial population size should be smaller 
than the carrying capacity (Hui 2006, Argasinski and Kozlowski 2008, Argasinski 
and Broom 2012) leading to 



where suppression (1 — n/K) describes newborns survival. This implies impor- 
tant distinction between newborn candidates introduced to the environment (de- 
scribed by b) and recruited newborns, survivors of density dependent stage (de- 
scribed by 6 ^1 — j^^)- This problem was emphasized by Kozlowski (1980) for the 
first time, but surprisingly this paper did not get as wide an appreciation as it de- 
served. Then it was reinvented by Ginzburg (1992), but ^ was rejected there as it 



(3) 
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"disagrees with our intuition about unchanging equilibrium " . Hui (2006) argued, 
against Ginsburg's claim, that is the proper approach and should substitute 
for ([2])- The discussion started by Ginzburg also did not receive wide attention. 
Without knowledge of this discussion, ([3|) was applied to avoid the suppression of 
selection after reaching the equilibrium size caused by (Argasinski and Kozlowski 
2008) and is a cornerstone of the ecologically realistic approach to dynamic evo- 
lutionary games (Argasinski and Broom 2012). Then ([3]) was mentioned as an ex- 
ample of the proper mechanistic approach at (Geritz and Kisdi 2012), but not as 
the general alternative to ([2]). However, we believe that ([3]) deserves much stronger 
attention from a general audience. Although ^ has been applied in complex selec- 
tion models (Argasinski and Kozlowski 2008, Zhang and Hui 2011, Argasinski and 
Broom 2012), the selection consequences of this approach were not rigorously ana- 
lyzed, since previous papers (Kuno 1991, Ginzburg 1992, Gabriel 2005, Hui 2006) 
were focused on population density dynamics and ecological aspects. The selection 
mechanism induced by ([3]) is very interesting and will be analyzed in the further 
sections. 

3. The population in equilibrium 

We can calculate equilibrium size, by setting the right hand side of Equation ([3]) 
to be equal to 0, which gives either n = or 



Note that for positive n, the condition b > d should be satisfied. After substitu- 
tion of n into the logistic coefficient (1 — n/ K), we obtain the equilibrium newborn 



(4) 




survival d/h. This is reasonable; due to the turnover of individuals we have nb new- 
borns and nd dead individuals per unit time. Thus nd/nb describes the number 
of newborns competing for single nest site vacated by a dead individual. Only one 
newborn can settle in a single place, thus each newborn can survive with probabil- 
ity d/b. This newborn survival should be valid for any density dependent mortality 
acting on juveniles, not only for logistic suppression, because only in this case does 
fertility equal overall mortality. 

4. The case of multiple individual strategies 

Assume that there are different individual phenotypes i = 1, . . . ,H each charac- 
terized by reproduction bi and mortality dj. Thus every strategy is described by a 
two dimensional vector Vi = [bi, di] G ([0, oo) x [0, 1]). 

Denoting n = J2i ^^'^ Qi = —,we can describe the following dynamics: 



n 



(5, f .„.(,(,-ii)-..). 



Thus 



i i \ i i / 



givmg 



where b{q) = Yliili^i ^"^^ = Yliili'^i- We can easily calculate that in this 
case, instead of reaching the stable equilibrium, population size converges to the 
stationary density manifold (Cressman et al 2001,Cressman and Garay 2003a and 
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b) 



(7) 



which form is conditional on the strategy frequencies. Thus we introduced the 
diversity among individual strategies to our model. Newly produced offspring of 
the carriers of the different strategies forms a pool of candidates from which a 
randomly drawn individuals will be recruited to settle in the available nest sites. 
This is the core of the " nest site lottery" mechanism which will be analyzed in the 
next sections. 



Here we will introduce an important characterization of the population dynam- 
ics. We shall define the function L{v) = h/d for a single strategy v. L describes 
the number of newborns per single dead individual, which we shall refer to as the 
turnover coefficient (for the relationship of the turnover coefficient with lifetime 
reproduction, see Discussion). Surprisingly, a similar coefficient describing the en- 
ergy allocation into reproduction divided by mortality can be found in life history 
papers (Taylor and Williams 1984, Kozlowski 1992, Kozlowski 1996, Werner and 
Anholt 1993, Perrin and Sibly 1993, for overview Kozlowski 2006). Analogously, 
for a mixture of strategies where v{q) = gjU, = [6, d\ is the average strategy con- 
tained in the convex hull of the strategies Vi (see figure 1), we define L{v{q)) = b/d. 
Thus 



5. The turnover coefficient L 
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which is a weighted average of the LiS and yi — qidi/Y^ qjdj describes the fraction 
of i strategists among individuals dying during a small time interval At (according 
to Appendix A, IA.2p . This is thus the average Li among dead adult individuals. 
The L-function can be useful in describing the multiplicative newborn survival 
(recruitment probability) because after substitution of the stable population size n 
into the logistic suppression coefficient we obtain: 



which can be denoted as 1/L{v{q)). 

If there is any variation in the i^s, then we have that d/b lies strictly between 
the smallest and largest values of L, Ln^n < b/d < imax- L{v) describes the 
number of newborn candidates produced per single dead individual for the strategy 
V = [b, d], during At. When the strategic argument is the averaged vector describing 
a population with a mixture of strategies, then the value of L is the average number 
of newborn candidates produced per single individual in this population. When 
the population is in size equilibrium (at the stationary density manifold), then the 
newborn survival component can be described by the value of L of the average 
population strategy, thus it becomes frequency dependent. 



Now we can describe the selection process realized by "nest site lottery" mech- 
anism. Thus we should describe population in terms of the strategy frequencies 
Qi and the population size n In our model, in the general case, newborn survival 
^1 — — ^ is the phenomenological function, linear with respect to the fraction of 
free nest sites. However, at the stable size manifold it has the form We can 



(9) 




6. Selection dynamics 
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assume that the strategies are enough close to each other that a separation of 
timescales between fast n dynamics and q dynamics occurs. Then we can assume 
that the selection occurs on the stationary size manifold. Thus using ([9]), we can 
write the selection dynamics from ([5]) as 

(10) —nr='^i\°i 



dt 'V Hq) 

Because the average growth rate on the stable size manifold is zero then the 
equations (jlOl) can be substituted by the replicator dynamics (see Appendix A) 

Therefore the growth rate of the j-th strategy becomes a function of the strategy 
frequencies q (frequency dependent).' 

and by equation ([6]) the population size manifold is described by 



biq). 

Note that the growth rate function M describes a mixture of all mortality and 
fecundity components, not only the density independent mortality d and fecundity b 
as in the Malthusian parameter r. The growth rate is positive when L{vi) > L(v{q)), 
which implies that 



Thus there is a threshold between regions of growth (strategies with reproductive 
surplus) and reduction (strategies with death rate exceeding birth rate) which has 
the linear form bi — L{J2i <liVi)di (see Figure 1). The threshold describes the set of 



12 



strategies for which the growth rate M{vi, qi) equals 0. Frequencies qi of strategies 
Vi with a greater value of L than the average strategy ^ ■ qiVi will increase under the 
replicator dynamics. In effect the averaged strategy shifts towards those strategies 
because it is a linear combination of the strategies present in the population. This 
implies an increase of L of the average strategy (see Figure 1). However, among 
growing strategies, the greatest growth rate is by the strategy with the greatest 
Malthusian parameter M{vi,qi). Frequency dependence induces an increase of the 
slope of the threshold which eventually leads to the selection of the strategy with 
the greatest L, which confirms the result of Mylius and Diekmann (1995). Note 
that their second result, that density dependent adult mortality leads simply to r 
maximization as in unlimited growth models, directly comes from the independence 
of replicator dynamics from background fitness. 



We have seen that evolution leads to the fixation of the strategy with the largest 
value of L, Lmax- What if there is more than one such strategy? The following 
question arises: is there selection between strategies with the same L? We can 
show this by applying a multipopulation game-theoretic approach (Appendix A 
and Argasinski 2006) and divide strategies present in the population among sub- 
populations with the same L, but diff'erent ds (Appendix B). In effect we obtain 
the following equations 



7. Multiple strategies with Li = L, 



(12) 




(13) 
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describing the selection of strategies within the j-th L-class (equation (fT2|) ) and 
related frequencies between L-classes described by gj. Thus selection between 
L-classes is driven by the first bracketed term from Equation (fT2|) and affects both 
intra- and inter-group dynamics. However, there is selection inside each L-class 
toward greater d. When suboptimal i-classes are outcompeted, the intrinsic selec- 
tion driven by the bracket (^d{ — XI ^uj'^iu) ^'l^o suppressed. The form of Equation 
([T3l) shows that among growing L-classes, those with smaller L can grow faster than 
those with larger, due to a greater J^lw^-L^ until they fall under the L-selection 
threshold. 

We are interested in the question which strategy (strategies) will dominate the 
population in the long term. In particular, qi is always positive if strategy i has 
Li = imax, whenever there is variation in the L values in the population. Thus 
the proportions of such strategies increase; but also, following Equation (IT^ . the 
strategies out of these with the largest values of di increase the fastest. Thus if 
there is either repeated small mutations involving strategies with Li < Lmax or a 
constant low level of mutation involving a mix of strategies making L{v{q)) < L^i^^, 
the population will evolve to the strategy out of those with Li — Lmax such that 
di takes the largest value. Thus repeated mutations or invasions of suboptimal 
strategies induce selection towards maximal d among Lmax strategists. It is easy 
to show that in the absence of density dependent suppression, this strategy has the 
greatest r but only among Lmax strategies, since bi — diL^^^- The strategies from 
other L-classes can have even greater r, but they will be outcompeted by mechanism 
described by bracketed term in equation . In the case when population consists 
only of the Lmax individuals, the same outcome can be caused be repeated ecological 
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catastrophes leading to the decrease of the population size. Then the strategy with 
the greatest di will have the greatest growth rate during the growth phase of the 
population. 

However, let us focus on the evolution of the system under the replicator dy- 
namics in a single particular "turn", during which no mutation occurs. Suppose 
that there are precisely / strategies in the imax-class. Assume that initial state of 
the Lmax-class is described by vector q = (g™'''^(0), . . . , c^'^^i^)) and initial relative 

size 5'max* 

Theorem 1 

The replicator dynamics converges to the vector q = (g™'''^(0)5max(0)A'*^ , 
where A is a constant that satisfies the equation 5(max(0) 9j(0)A'^* = 1 
For the proof see Appendix C. 

The above theorem shows that the restpoint describing frequencies among Lmax 
strategists is fully determined by the initial state of the Lmax-class and its relative 
size. Despite frequency dependence, independently of the initial frequencies of the 
other strategies (note that in our model there are no direct interactions between 
individuals). Thus calculation of q reduces to finding of the apropriate value of A. 
Note that the rest point q can be interpreted as the state of the whole population 
in general coordinates and the final stateof the Lmax-class (then flmax = !)• If 
initial 5max(0) = 1 then obviously A = 1, thus parameter A can be described as 
the inflation coefficient because it inflates the frequencies to sum them to one and 
compensate the impact of 5max(0). 



...,gr"(0)5ma.(0)A''O, 
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8. The monomorphic resident-mutant case 

We can simplify the above model by assuming a monomorphic population in- 
vaded by a rare mutant; thus this resembles the classical ESS approach (Maynard 
Smith 1982, in the context of life history evolution Charlesworth & Leon 1976, 
Mylius and Diekmann 1995). In the limit as the change of the strategy trait tends 
to zero, we approach the method known as invasion analysis which is the cornerstone 
of adaptive dynamics (Dieckmann and Law 1996, Metz et al. 1996, Geritz et al. 
1998, Dercole and Rinaldi 2008). Using Equation ([TT|) . the resident Malthusian pa- 
rameter is zero and the rare mutant Malthusian parameter is {bmut/bres)dres — dmut 
which must be positive to invade the population. Thus and the equilibrium popu- 
lation size increases. 

To be an ESS itself, the "mutant" population should be stable against the pre- 
vious resident, and thus: 



I I ^ n ^ ^res ^miit 



mut (J'res 



which is the same condition. 

Therefore we have obtained for the "nest site lottery" mechanism, the classical 
result that under limited growth only lifetime reproduction is maximized and there 
is no selection pressure on the lifespan. However this occurs only in a monomorphic 
resident-mutant model. 
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9. Discussion 

9.1. How the nest site lottery works? We started from the basic population 
growth equation which is the cornerstone of the framework underlying evolution- 
ary game theory and replicator dynamics (Maynard Smith 1982, Cressman 1992, 
Hofbauer and Sigmund 1988 and 1998) and its more ecologically realistic exten- 
sions (Cressman and Garay 2004, Argasinski 2006, Argasinski and Kozlowski 2008, 
Argasinski and Broom 2012). We presented an analysis of the dynamics of the 
mechanism inducing frequency dependent selection toward the strategy maximiz- 
ing the turnover coefficient L{vi). 

This phenomenon can be explained mechanistically. All newborns introduced 
to the population in the same time moment form a pool of candidates. Each 
newborn is independently drawn to survive (find a nest place), thus the strategy 
maximizing the number of newborns (trials) maximizes the fraction in the pool of 
candidates and in effect the amount of survivors. However, every dead adult can 
be substituted by any other strategy, thus each death is an additional free place in 
the lottery. Thus it is profitable for the strategy carried by some subpopulation to 
maximize the number of trials (newborns) per single offered place (dead adult). In 
addition, we have shown that among strategies with the largest value of the turnover 
coefficient Lmax there is a selection pressure toward the strategy with the greatest 
d. This is intuitive from equation ()lip . because for strategies with the maximal 
number of newborn candidates produced per each dead adult (i.e. maximizing the 
bracketed term in equation (|11|) ). the growth rate will increase with the number 
of dead adults (described by the fraction qidi in the equations pT]) ) since each of 
them will be exchanged for Lmax newborns in the pool of candidates. Note that 
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this provides a gene centered mechanistic explanation of the phenomenon which 
can be naively interpreted in terms of group selection and altruistic "sacrifice" of 
adults, to release the nest sites for juveniles. We show that it is the tradeoff between 
maximization of the number of candidates in the nest site lottery and survival of 
parental individual. 

9.2. Importance of the nest site lottery mechanism. The model presented 
in this paper is as simple as possible, to emphasize the mechanistic aspects of the 
analyzed phenomenon. For example, there are no direct interactions between in- 
dividuals as in game-theoretic models. However, it was shown that the impact of 
density dependent factors (thus also the mechanism described in this paper) can 
significantly affect and alter the outcomes of game-theoretic models (Argasinski 
and Kozlowski 2008, Argasinski and Broom 2012). This is caused by the feedback 
driven by the fact that on stable size manifold every newborn should find a new 
nest site vacated by a dead adult. Our model is the simple case example of the 
one- dimensional monotone density dependence acting on effective birth rate (more 
on this and other general cases can be found at Metz et al 2008). Our model has 
extremely simplified age structure consisting only of juveniles and adults (more on 
the limitations of pure age-dependent models can be found at Metz and Diekmann 
1986). However, it is possible that other environmental feedback loops of the same 
type than in our model may induce similar selection mechanisms. Maybe phenom- 
ena similar to those revealed by our simple model, can be found in other, more 
general or different specific models related to general population dynamics, life his- 
tory evolution, adaptive dynamics or population genetics. This can be subject of 
the future research. 
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The mechanism shown in this paper supports the intuition underlying the r and 
K selection theory (McArthur and Wilson 1967), that natural selection favours 
different strategies in growing populations than in populations with suppressed 
growth. The theoretical and methodological aspects of this approach were criticized 
(Barbault 1987, Getz 1993, Stearns 1977), however as an intuition it still seems to be 
relevant (for modern approaches see for example Metz et al. 2008). The alternative 
for r and K approach is the life history theory (Roff 1992, Stearns 1992), where 
the problem of different selection mechanisms in limited and unlimited populations 
also exists. 

9.3. What is maximized by natural selection, and when? The exact mean- 
ing of "fitness" is a subject of endless discussion (Metz et al. 1992, Kozlowski 1993, 
Mylius and Diekmann 1995, Brommer 2000, Dieckmann and Metz 2006, Metz et 
al. 2008, Roff 2008). Basically "fitness" can be defined as the instantenous growth 
rate or invasion exponent (Metz 2008c). However, if eco-evolutionary feedback is 
of particulary simple kind, the optimization approach can be applied (Metz et al 
2008, Gyllenberg et al 2011) where some "fitness measures" or "proxies" are maxi- 
mized. There is a widely known fact in life history theory that in a population with 
unlimited growth, the Malthusian parameter r is a proper fitness measure, while 
on stable size manifold, lifetime production of newborns (before juvenile mortality 
selection) Rq is the correct measure. In (Mylius and Diekmann 1995) there is a 
statement that the invasion fitness method (Metz et al. 1992) suggests that Rq and 
r are necessarily both maximal at the ESS (although this statement is unclear since 
some strategies can maximize Rq and others can maximize r). Our results support 
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this claim, and show that an analogous mechanism can act in population dynamic 
and game-theoretic models. 

Here an important claim is that of Brommer and Kokko (2002), who say that Rq 
is a rate independent reproductive measure which does not depend on the timing of 
reproductive events. This is because Rq is described on the lifespan timescale, not 
the population dynamic timescale like r. Despite simplified form our model can be 
useful illustrative example for this problem. How can we use a lifespan perspective 
in our approach? At first, assume that the population growth is unlimited (every 
newborn candidate can find a nest site). Then demographic parameters b and d are 
constant. When we consider only the death component of the Malthusian equa- 
tion then we obtain the equation h = —nd; thus this is an exponential decay with 
decay constant d, and so the average lifetime of the individual is r = 1/d. Then 
r = {L — 1) /t and L{v) is the lifetime reproduction Rq (or the average Rq among 
individuals dying during At in a population described by v). Therefore, the for- 
mula r = {L — 1) /t shows how the growth rate is affected by lifetime reproduction 
in the case of a non-age structured population with unlimited growth. It shows 
that an individual should basically replace itself, but for strategy growth rate to be 
positive requires a reproductive surplus during a lifetime. Now introduce the limi- 
tation of the nest sites causing density dependent selection. Our result shows that 
under density dependence, the growth rates of the strategies are affected by the 
frequency dependent surplus reducing mechanism, described by newborn survival 
1/L{v) which is the function of the other strategies present in the population. In 
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effect r{v) = {L{v) / L{v) — 1) /t and i?o = L(v)/L{v). Thus the strategy maximiz- 
ing the Ufetime production of newborn candidates L{v) will maximize Rq and the 
bracketed term of r{v). 

9.4. Conclusion. Our simple model suggests an insight into the mechanistic na- 
ture of this phenomenon and has serious interpretational consequences. It clearly 
shows that this problem should not be formulated as the alternative: evolution 
maximizes r OR Rq. In our simple model, when the population reaches a stable 
size manifold, then a mechanism that modifies the r's of competing strategics, which 
are no longer constants, emerges to select the strategy with maximal Ro, or with 
maximal r, among multiple strategies with maximal Rq. Thus, our model suggests 
the existence of another fitness measure which is the combination of i?o and r, if 
our reasoning holds in age structured and other more complex models. However, 
it will be not a fuction which should be maximized, but a two staged procedure. 
First stage should identify the strategies maximizing the turover coefficient, while 
second stage should find strategies with the greatest r from strategies chosen in the 
first stage. We note that our analysis is a simplification, and whereas Rq = L{vi) in 
the models without age structure as presented in this paper, this is not necessarily 
satisfied in age structured models. This should be the subject of future research. 
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Appendix A Multipopulation replicator dynamics. 

Assume that we have u individual strategies. Standard replicator dynamics can 
be derived by rescaling the Malthusian growth equations ^ = UiMi to the related 
frequencies Qi = rii/ J2j '^j which leads to equations ^ = Qi [Mi — M] (where M = 
J^jQjMj)- Those equations describe the evolution of strategy frequencies in the 
unstructured population. However, we might be interested in the modelling of the 
structured population divided into subpopulations such as different sexes, species 
etc. Assume that we want to decompose an entire population into z subgroups. 
Define 

(A.l) = \ki,...,ki^] 

as a vector of indices of strategies exhibited by individuals from the i-th subgroup 
(fcj G {!,...,«}, and uj is the number of strategies in the i-th subgroup). For 
example the notation k'^ = [1,3,5] means that, in the second subgroup there are 
(only) individuals with strategies 1,3 and 5. Every strategy should belong to a 
single unique subgroup. Then according to Argasinski (2006), by the following 
change of coordinates 



(A.2) q' = [ql...,ql] 



we obtain a distribution of relative frequencies of strategies in the i-th subpop- 
ulation. The distribution of proportions between subpopulations has the form 



(A.3) 



9 = 
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where gj is the proportion of the j-th subpopulation. Every decomposition into 
subpopulations can be reduced again to a single population model by the opposite 
change of coordinates q{g, q^, q^) where 



When we apply the above transformations to the replicator equations, we obtain 
a set of equations that describes the dynamics inside the subpopulations (intraspe- 
cific dynamics). When the set of strategies in each subpopulation is characterized 
by the vector of indices , then the system of replicator equations will be: 

i = l,...,Uj - 1, j = l,...,z 



(A.5) 



aq; 
dt 



--Qi 



(A.6) ^=^4M«_M] s = l,...,z-l 

where = it^i is the mean fitness in the s-th subpopulation and M = 
9sM^. In practical applications of this method to the modelling of biological 
problems, the replicator equations can be defined on the decomposed population. 
This will simplify the formulation of the model, because when strategies are initially 
assigned to subpopulations, there is then no need to change their indices. The choice 
of subpopulations is arbitrary and depends on the biological assumptions underlying 
the analyzed problem. For example, the entire population may be divided into 
two competing subpopulations of hosts and parasites or prey and predators. On 
the other hand, it may be divided into two subpopulations of males and females, 
when interspecific dynamics will describe the evolution of the secondary sex ratio. 
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and intraspecific dynamics will describe changes of the frequencies of strategies 
inside the male and female subpopulations. The subpopulations can be divided into 
subsubpopulations, and the entire population may be transformed into a complex 
multilevel cluster structure. However, all these structures are equivalent to a single 
population replicator dynamics model. 



Appendix B Derivation of equations (|12|) and (|13|) describing selection 
strategies inside L-classes and change of sizes of i-classes 

Let us assume than we have m such classes. Then the initial system of the repli- 
cator equations can be transformed into two sets of differential equations. Firstly, 
the within L-class dynamics (according to IA.5|) : 

(B.l) ^^ql^Mivi)~M^), 

where qf is the proportion of the z-th strategy in the j-th L-class and = 
'Ewli^'HO = 'Ewli'^i {L{vl) I L{v{q)) - 1). Secondly, the between L-class dy- 
namics (according to IA.6P : 

(B.2) ^ = [W - M) , 

where is the proportion of j-th L-class and Af = 0, since the population is on 
the stable size manifold. Since for all strategies (for all i) from the same L-class 
the term {L{vl) / L{v{q)) ~ 1) in B.l is the same within the class, the lower index 
i (representing the strategy number within that class) can be neglected. After 
substitution of the respective formulae into Equations (|B.1[) and ()B.2p , we obtain 
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the equations ([T2|) and (|T3l): 



dt \L{v{q)) 



Appendix C: Proof of Theorem 1 

From we have that 



% _ , f L{Pi) _ . 



and so 



^ \L(v{q)) J diQi dt 

Consider any pair of strategics Vi — [fe^, d^ and Vj = [6j, dj\ from the same L-class 

(i.e. L (vi) = L {vj)). Using (jC.l|) we obtain 

1 dqi 1 dqj 
diqi dt djqj dt 

1 f 1 
-dqi = / 3 dqj + C 



djqj J djqj 



(C.2) In^^ln^^^^ 

di dj 



Considering i = in equation (|C.2p we obtain 

lng,(0) Ingj(O) 



(C.3) C = 



di dj 



Combining (IC.3|) with (jC.2[) we obtain 

lng,(t) - lng,;(0) _ Ingj (t) - Ingj (O) 



\<lo{t)J \qoit) 
Equation (|C.4|) holds for any pair i,j so that 

for some A(t). Letting A = linif^oo A(i) gives 
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(C.5) (7, = hm q,{t) = q,{0)X 



di 



We know that J2i Qi — li thus for at least one L-class respective X{t) should not 
converge to 0. System p2|) . (fT3|) shows that it will be imax-class. However, above 
reasoning was described in the coordinates describing the strategy frequencies in 
the whole population (a metasimplex coordinates, Argasinski 2006). According to 
(IA.4I) . qu{0) can be described in the coordinates of the system ([T2| and (fT3|) and 
after change of the indexes u — k\ where I is the index of the L-class and a is the 
index of the strategy within this L-class, we have qi^i — giq\- Thus the restpoint 
will contain only the L-maximizing strategies i.e state of the Lmax-class will be 
equivalent to the state of the whole population (i.e according to \KA\ .a^a^ = 1 
and gfcmax — q™'^^), but frequencies gfc>"'>x(0) will not sum to 1. However, from (jA.4p 
we have gfcmax(O) — 5max(0)g™^''(0). Then (jC.SP for the imax-class can be presented 



(14) gr'^ = gr'^(0)5.„ax(0)A^". 
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Figure legends: 
Figure 1 

The figure shows the plot of the i-surface over the strategy space. The indi- 
vidual strategies and i-classes are shown on the strategy space. The L-class of 
the populations average strategy is the border between regions of growth and de- 
cline. The growth of the frequencies of the strategies from the growth area induce 
counterclockwise movement of the i-class of the average strategy, leading to the 
selection of L-maximizing strategies. However, among L-maximizers, selection for 
greater mortality d can be induced by the introduction of sub-optimal mutants. In 
this case the strategy lying on the maximal L-class line, farthest from zero, wins. 



